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Abstract 

We apply (Fractional) Analytic Perturbation Theory (FAPT) to the QCD analysis of the non¬ 
singlet nucleon structure function F 2 (x,Q 2 ) in deep inelastic scattering up to the next leading 
order and compare the results with ones obtained within the standard perturbation QCD. Based 
on a popular parameterization of the corresponding parton distribution we perform the analysis 
within the Jacobi Polynomial formalism and under the control of the numerical inverse Mellin 
transform. To reveal the main features of the FAPT two-loop approach, we consider a wide range 
of momentum transfer from high Q 2 ~ 100 GeV 2 to low Q 2 ~ 0.3 GeV 2 where the approach still 
works. 
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I. Introduction 


QCD analysis of deep-inelastic scattering (DIS) data provides one with new knowledge of 
hadron physics and serves as a test of reliability of onr theoretical understanding of the hard 
scattering of leptons and hadrons. At large momentum transfer q, — q 2 = Q 2 1 GeV 2 we 
have the reliable description of DIS that is based on the twist expansion and “factorization 1 ' 
theorems. At small (moderate) transfer Q 2 < 1 (a few) GeV 2 this QCD description faces 
two main problems: (i) the high twist corrections to the leading twist contribution become 
important but remains poorly known; (ii) perturbative QCD (pQCD) becomes unreliable 
due to the fact that the QCD running coupling a s (Q 2 ) grows and “feels” infra-red Landau 
singularity appearing at the scale Q ~ Aqcd ~ of a few tenth of GeV. We discuss in this 
paper a solution of the last problem by applying to DIS analysis a nonpower perturbative 
theory whose couplings have no singularity at Q 2 > 0 and whose corresponding series possess 
a better convergence at low Q 2 . 

A widely used approach to resolve the aforementioned problem is to apply the Analytic 
Perturbation Theory (APT) developed by Shirkov, Solovtsov et al. [l!-@]. There, the running 
QCD coupling a s (Q 2 ) = a s (Q 2 )/4:7T of pQCD is transformed into an analytic (holomorphic) 
function of Q 2 , al(Q 2 ) H > Ai(Q 2 ), APT coupling. This was achieved by keeping in the 
dispersion relation the spectral density p[ pt ^(<7) = Im a s (Q 2 = —a — ie)/ tt unchanged on the 
entire negative axis in the complex Q 2 -plane (i.e., for cr > 0), and setting it equal to zero 
along the unphysical cut 0 < Q 2 < A 2 . In the framework of APT the images A n (Q 2 ) of 
integer powers of the originals a™(Q 2 ), a™ A n following the same dispersion relations were 
also constructed. At low Q 2 the couplings A n (Q 2 ) change slowly with Q 2 in contrast with the 
original a™(Q 2 ) behaviour while at high Q 2 A n (Q 2 ) —> a™(Q 2 ). Later, the correspondence 
a v s i — y A u was extended to noninteger powers/indices v in [6-lo( and was called Fractional 
APT (FAPT), which provides the basis for application to DIS. In this respect let us mention 


a recent papers [11] where the processing of the DIS data has been performed in FAPT in 


the one-loop approximation and the reasonable results for hadron characteristics has been 
obtained. 

Various analytic QCD models can be constructed, and have been proposed in the liter¬ 
ature, among them in Refs. 13—1191]. These models fulfill certain additional constraints at 



For further literature on various analytic QCD models, we refer 
Some newer constructions of analytic models in QCD of 


low and/or at high 
to review articles [1 

Ai(Q 2 ) include those based on specific classes of (5 functions with nonperturbative contribu- 

and those based on modifications of the the 
ze)/7r] at low (positive) cr where pi is 


-cr 


tions [22] or without such contributions [23H25 
spectral density p^ i->- p, [pi[cr] = Im A±(Q 2 = 

parameterized in a specific manner by adding two positive delta functions to p ( j pt ' ) , cf. 

The possibility to extend the DIS analysis formally in the whole Q 2 range together with 
the effect of slowing-down of the FAPT evolution of the parton distribution functions (PDF) 
in the low Q 2 region are attractive phenomenological features of FAPT. A number of works 
deal with this task in a naive form [27], where the authors show that at very low Q 2 and 
Bjorken variable x APT agrees with experimental data. Besides, the applicability of the APT 
approach was analyzed in the Bjorken polarized sum rule [28| confirming that the range of 
validity of APT is down to Q ~ Aqcd — 350 MeV, as compared to experimental data. The 
common feature of these works was taking into consideration some nonperturbative effects 
against the background of APT, i.e., higher twists in [27l-29| or an effective constant gluon 
mass in 


30 . 
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The basis for applying FAPT to low energies in this approach is the factorization theorem 
that allows one to shift the frontier between the perturbative and nonpertnrbative effects 
via the variation of the factorization scale. Therefore, we shift the range where perturbation 
series is applicable in FAPT, as it was demonstrated in [28] (see reviews of this issue in [10], 
where this phenomenon was also discussed for pion form factors). 

Our goal here is to elaborate a general scheme of DIS data processing in the framework 
of FAPT taking as a pattern the DIS analysis at NLO. In this respect the discussion here 
can be considered as an extension of the partonic results of the article [nj on the higher- 
loop level. We shall focus on the specifics of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution of the PDF /(x; p 2 ) in FAPT. We involve into consideration the co¬ 
efficient function C(x,a s ) of the process and compare the final result with a similar one 
in pQCD. An important problem of higher twist contribution remains untouched here, but 
higher twist effects can be taken as an unknown function h(x), i.e. h(x)/(F Ilf, or as a 
constant jUj] p 4 /Q 2 , or as an effective sum of all twists contributions in [29|. We stress 
that HT effects are only indirectly affected by the analytization procedure. The behavior 
of HT will be given by the fit of experimental data together with the corresponding parton 
distribution functions fill . 27 . 28]. Besides, in [28[ the authors included more terms in the 


HT expansion and demonstrated that they are essentially smaller and quickly decreasing. 
Because of this (theoretically) unknown behavior we avoid this problem since we pretend 
to provide perturbation tools how to deal with FAPT, while the pure phenomenological 
analysis is transferred to future investigations. 

Let us recall that the DIS analysis can be performed in a few different ways: one of them 
is provided by the Mellin moments defined via inelastic structure functions (SFs) F(x,Q 2 ), 


M(n,Q)= / dxx 11 l F(x,Q ), (n — 1,2,3,...). 


(i.i) 


The second approach is based on the direct application of the DGLAP integro-differential 
evolution equations [31] to PDF /, while the observable SF is the Mellin convolution of 
the coefficient function and PDF, F = C * f. The third approach makes use of the Jacobi 


Polynomial expansion method 32]. Just this method will be used in this work. 


The paper is organized as follows: in Sec. HU we present a theoretical background where we 
describe the Jacobi Polynomial (JP) method and how to calculate free parameters in order to 
obtain the nonsinglet structure function. In Sec. IIIII we briefly describe the FAPT approach 
and derive the DGLAP evolution for the moments M{n , Q 2 ) in FAPT. We present in Sec. IIVI 
the free parameters obtained in the analysis of the so called MSTW parameterization, see [33| 
for details, and the nonsinglet SFs at LO. Section HVl contains the results of the analysis of 
numerical realization of the FAPT evolution and the comparison with the results of analogous 
calculations in pQCD. Finally, in Sec.|V]we summarize our conclusions. Important technical 
details including new findings are collected in four appendices. 


II. Jacobi Polynomial expansion for DIS analysis 

We shall focus here on nonsinglet (NS) structure functions, Fns(x,Q 2 ), with their cor¬ 
responding Mellin moments Mxs(n, Q 2 ) (via Eq. (11.ip ) to avoid technical complications of 
the coupled system solution in the singlet case. The PDFs / P (x,p 2 ) are universal process- 
independent densities explaining how the whole hadron momentum P is partitioned in x- P, 
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i.e., the momentum carried by the struck parton (see, for instance j34|). The ^-dependence 
of PDF is formed at a hadron scale of an order of P 2 by nonperturbative forces, while its 
dependence on factorization/renormalization scale p 2 can be obtained within perturbation 
theory. 

A brief description of the evolution of the Mellin moments in pQCD, up to NLO, is out¬ 
lined in Appendix [B] as well as the theoretical background with our notation and conventions. 
We consider the scale Ql as a reference scale for the solution of the evolution equation (IBGaj) 
where the PDFs are regarded as functions of x and the parameters are fixed by comparison 
with DIS data. In particular, we use here the data-based MSTW PDFs (see [33j, where 
Ql = 1 GeV 2 ). Namely, 

xu v (x,Ql) = A u x m (l - x) m (l + e u \fx + 7 „x), (2.1) 

xd v (x, Ql) = A d x m (l-x)' u (l + e d y/x-\-'y d x), (2.2) 

where the values of A ud , g k (k = 1,..., 4), e U:d and 7 , t/ / can be found in j33|. We use only the 
valence quark PDFs because the NS PDF / N g can be expressed as fws( x j Q 2 ) = w„(x, Q 2 ) — 
d v (x,Q 2 ) (see Appendix iBl for details). The NS SF F 2 (x) — (C * /ns) (x) is represented as 
the Mellin convolution of coefficient function C of the process and the corresponding PDF 
/ns- The F 2 (x) can be expanded in the Jacobi Polynomials Q%f(x), which was developed in 
Refs. [32], in truncating the expansion at n — N max , where the method converges (see, for 
review [35[): 

Nmax 

F 2 (x,Q 2 -N max )=uj a ^x) Q^(x)J2^ n \a,/3)M NS (j + l,Q 2 ). (2.3) 

n =0 j =0 

Here M^(n, Q 2 ) are the Mellin moments of nonsinglet SF calculated explicitly in Eq. (1B6alk 
uj a/3 (x) = x a (l—xY is the weight function and the parameters a, fd will be obtained by fitting 
to the data. The Jacobi Polynomials Q^(x) are defined as an expansion series by means of 

k 

efW = E c 'f , ( a ’^' < 2 - 4 > 

3=0 

They satisfy the orthogonality relation 

f u a P{x)Of{x)Qf(x) = 8 kl . (2.5) 

Jo 

Another way to obtain SF F 2 (x,Q 2 ) is to take the inverse Mellin transform M.~ l under 
the moments Mns(^, Q 2 ) (he., the inverse of Eq. (Il.ip ). Choosing a convenient path of 
integration one obtains for F 2 

-t /»c+zoo 

F 2 (x, Q 2 ) = M^ 1 {M NS (n,Q 2 )j = — / x~ n M NS (n,Q 2 )dn, (2.6) 

Z7TZ J c—ioc 

here we take the path along a vertical line Re(n) = c. We perform the “exact” numerical 
Inverse Mellin Transform, further comparing the results with the Jacobi Polynomial method, 
only at the one-loop level due to technical limitations. In this way, we estimate the accuracy 
of the applied polynomial method, the results of this numerical verification are outlined in 
Appendix O 

Let us finally mention that one can take SF F 2j instead of the NS F 2 to consider, e.g. , the 
neutrino DIS results of the CCFR collaboration, like it was started in [ll|. This replacement 
will lead to only minor changes of technical details in the procedure elaborated below. 
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III. FRACTIONAL ANALYTIC PERTURBATION THEORY and DIS 

It is known that the perturbative QCD coupling suffers from unphysical (Landau) sin¬ 
gularities at Q 2 ~ A 2 . This prevents the application of perturbative QCD in the low- 
momentum spacelike regime and, in part, impedes the investigation of high twists in DIS. 
Our goal here is not to discuss the motivation and complete construction of FAPT, which 
couplings A u are free of the aforementioned problems, but present to reader illustrations 
of the properties of this nonpower perturbation theory that are important for further DIS 
analysis. 


A. Elements of FAPT 


Application of the Cauchy theorem to the running coupling a"(Q 2 ) = {a s (Q 2 ) / A.n) 


established in and developed in joj fiol 


gives us the following dispersion relation 

Kallen-Lehmann spectral representation) for the images A ( ,P in the spacelike domain 


,or 




p { i\L a ) 

1 + exp(L - L a ) 


dL a , 


(3.1) 


(where L a = ln(cr/A 2 )) that has no unphysical (Landau) singularities. For the timelike 
regime analogous coupling reads 



p®{L a )dL 0 


(3.2) 


Here, Au(L) is the FAPT image of the QCD coupling a v s ^(L) in the Euclidean (spacelike) 
domain with L = ln(Q 2 /A 2 ) and the label l denotes running in the l -loop approximation, 
whereas in the Minkowski (timelike) domain, we used in (13. 2 p L s = ln(s/A 2 ). It is convenient 
to use the following representation for the spectral densities pp : 


R m lL) 


-Im (o, (J )(L - in))" = 

a,(i)(L -mt)|; ippL) 


IT (%(£))■'’ 

= arg (a s g,(L - nr)). 


(3.3) 


From the definition (13.11) and Eq. (l3.3l) it follow that there is no standard algebra for the 
images A u , i.e. A U A^ A u+I1 that justifies the name nonpower perturbative theory. 

In the one-loop approximation, the ipm, Rm has the simplest form, i.e., 


(f(i )(L) = arccos 


y/T ?T 


7T^ 


R(1){L) = PoV L 2 + 7T 2 . 


(3.4) 


Substituting Eq. (I3.4[) in Eq. (13.31) for p\ 1 } and then the result p^\L a ) in Eq. (l3.1|) . one repro¬ 
duces at Q 2 = 0 the well-known expression for maximum value of A^(L), A^\l = — oo) 

0, 



dL a 

A) (A 2 +vr2) 


1 

Po 


>A?(L). 


(3.5) 
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At the two-loop level, they have a more complicated form. To be precise, one gets 


1 1 

^1 + W- 1 {z w (L)Y 


(3.6) 


and 


R(2){L) 

^(2 ){ L ) 


ci(n f ) |1 + W-i(z w (L + Z7r))|, 

Re (1 + W- 1 (z w (L + in))) 


arccos 


R(2)(L) 


(3.7) 


with W_i(z) being the appropriate branch of the Lambert function, Z\y(L) = 
—ch 1 (n/)e _1 ' i / Cl(n ^\ c k {iif) = /3 fc (n/)//d 0 (n/) fc+1 , where /3 k are the QCD /3-function co¬ 
efficients and rif is the number of active quarks, see the expressions in Appendix |A] For our 

purpose we use here only the two-loop couplings like aiL , A u ■ Extensions up to four-loops 
can be found in {36}. 

Now we implement this formalism with the help of numerical calculation with the main 
module Mathematica package FAPT.m of [36} (confirmed by a recent program in [[ 37 }). 
According to this and using the corresponding notation from [36] in the RHS of Eqs. (13.81 
EM, we have 


AP(L) 

a <?(£,) 


AcalBarZ[L, n /, v\ 


[AnY 

UcalBarZ[L, nf , u] 


(47T ) v 


(l — 1 -j- 4; rif — 3 -j- 6) 
(Z = 14-4; n s = 3-5-6) 


(3.8) 

(3.9) 


For the coupling in pQCD we obtain 


a s (i)(L = \n(Q 2 / A 2 )) 


aBarZ [Q 2 , n f , A] (/ = 1 ^ 4 ^ 


(3.10) 


The correspondence between the pQCD expansion and FAPT one is based on the linearity 
of the transforms in Eqs. (13. II) and (13.21) . see [5]. This can be illustrated for the simple case 
of a single scale quantity D(Q 2 ,/j, 2 r ), calculated within minimal subtraction renormalization 
schemes and taken at the renormalization scale /i 2 R = Q 2 . The expansions for D and for its 
image D t-)- V are written as 


pQCD: D(Q 2 ) = d 0 <(Q 2 ) +J2 d n< + "(Q 2 ) 

n 

FAPT: V(Q 2 ) = d 0 A u (Q 2 ) + A [n+v) {Q 2 ) (3.11) 


at the same coefficients di that are numbers at n 2 R = Q 2 - 


B. FAPT for DGLAP evolution in NLO approximation 

We start with the well-known solution of DGLAP equation for the nonsinglet PDF /ns in 
NLO approximation. This solution is combined with the corresponding coefficient function 
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C(x,a s ) - the parton cross-section taken at the parton momentum xP. This is presented in 
Appendix iBl in the form of Eq. (1B6a| for the moments Mns of the NS SF F 2 . 

Rewriting Eq. (IB6af in the approximate form, i.e., neglecting the 0(a 2 s ) terms in the 
two-loop evolution factor, one arrives at the commonly used relation 


A/ns(^, Q 2 ) 


■yw+( c nJ(»)+| p(»)) «y )+i (Q 2 ) 

4f Al) + (4s(«) + |p(n)j «»S ( ” >+1 W§) 


(3.12) 


The use of FAPT will change in this scheme the sense of expansion parameters a s in accor¬ 
dance with (13. lip . An analogous evolution relation for the analytic images of the moments 
M n g, M ns i ^ AJns, can be obtained from Eq. (13.121) by replacing the powers ( a s ) v with the 
FAPT couplings A u (with v being here an index rather than a power) [8|] and reads 


kA , ^ A dL{n)(Q 2 ) + (4s fa) + ^n)+l(Q 2 a , 

Af N s(n,Q ) = —-T-V— - M 


+ (^Ns( n ) + % p( n )) ^j s(n)+1 m 


A) 


I (2) 


NS 


Ml). (3.13) 


The implementation of the proposed calculation in the form of 


is quite direct. The 


FAPT evolution relation (13.13(1 for the moments is the main result of the Section. Further, 
we shall use code (13.8(1 from |36] to obtain A (2> (L), 


1 s( 2 )(L)j numerically. 

The last approximation was taken up to 0(^4.d NS +i) since the contribution of the next 
term in the FAPT expansion in Eq. (l3.13p is negligible in comparison with the previous one 
(as we demonstrate in Appendix 0). This analytic version of the moment evolution does not 
face any problems at low energies due to the boundedness of couplings and rapid convergence 
of the FAPT series. 

In the absence of a fit of experimental data for the FAPT model we propose a relation 
for the initial moments at Q^: 


/ns(w, Ql) — 


Ql 


a s (Qo) dNs(n) + (C'nsW + wA n )) a s (<3o) dNs(n)+1 

A^nsi 


(n,Qo) 


•^*ra(n)(3o) + (CnS^ 77 -) + foP( n )j “^dNsW+l^o) 


(3.14) 


where the moment of PDF (see Eq. (lB3l) ) in pQCD stands in the LHS, while the moment for 
PDF in FAPT stands in the RHS of the second equation. In other words, we take the same 
initial PDF as in pQCD from the MSTW data for these both cases (in [ll| the parameters 
were taken the same since the difference between them was negligible). We can use either 
the Jacobi Polynomial expansion or directly the inverse Mellin transform ( Appendix O) . 


IV. Results of numerical analysis 


The accuracy of the SF approximation by a finite number of Jacobi Polynomials (trun¬ 
cated at N max ) depends on the choice of the weight-function parameters. Therefore, we test 
the nonsinglet SF, given by the MSTW data, by searching for the minimum of ( Q 2 = Ql): 


Xa,p 


ri(theor),Nmax i A ex P) 

^2 !r 2 



(4.1) 
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where we have used Eqs. (EUi ctncl (jBGaj) at Q" — Qo' Thus, wg have T 2 (x, ^q) = 
F ,2 Xp ' l (x, Ql) and from Eq. (I2.3[) p^ heor ^ Nm ™ Q2) = F 2 Vm “ ;<: (x, Ql). Then, we deter¬ 
mine the values of a and j3 that provide the best fit to the data for different values of 
Nmax • At the one loop level we find: N max = 13, a = 0.05, and /3 = 3.03 for \ 2 & 10“ 9 , 
whereas for two loops we get (for even PDFs only): N max = 13, a = —0.8, and f3 = 2.99 
for x 2 ~ 10 -9 . To evolve nonsinglet moments, we need to fix the values of the QCD 
scale A 12 (ti f = 3) in the leading and next-to-leading order, taken in [33;] from the com¬ 
parison with data, where ai lloop \Ql = lGeV 2 ) = 0.682 =>- A^n/ = 3) = 0.359 GeV 
and otf loop \Ql = 1 GeV 2 ) = 0.491 => A 2 {n f = 3) = 0.402 GeV. In the case of FAPT. 
the scales Af APT (n/ = 3) must be taken into account very carefully. The authors of fllj | 
fixed the A value directly from the comparison with the data in the leading order (where 
Qo = 3 GeV 2 ) and obtained Af APT (n/ = 4) = 0.275 ± 0.039 GeV that corresponds to 
A p APT( n ^ = 3) = 0.333 ± 0.050 GeV. We can see that the perturbative and the analytic 
values of A are close to each other at least inside the margin of errors. For this reason, we 
will take Ai, 2 (% = 3) ~ A PAPT (nj = 3) for simplicity (recalling that an appropriate value 
should be taken from the analysis of the experimental data but this goes beyond the scope 
of this work). The couplings in pQCD and in FAPT were calculated with the Mathematica 
package developed by Bakulev and Khandramai in j36j] where the heavy flavour thresholds 
were taken into account. 

Taking into account the above estimates of the initial parameters, we substitute 
Eqs. (I3.12p and (I3.13P into Eq. (12.3p . and obtain the evolution of SFs up to NLO in pQCD 
or FAPT, respectively. We show the final results of the evolution in Figs. [1] fusing for DIS 
the character interval 0.3 < Q 2 < 100 GeV 2 . In Fig. |T] we fix x at two different values: 






FIG. 1: Nonsinglet SF F 2 (x, Q 2 ) vs Q 2 at (a) LO and (b) NLO. The Bjorken x = 0.25 for (a.l) and (b.l), 
and x = 0.7 for (a.2) or (b.2). The solid line represents the FAPT results and the dashed line - the pQCD 
ones. 


x = 0.25 in (a.l), (b.l), and x = 0.7 in (a.2), (b.2), where (a) and (b) represent the LO and 
NLO results, respectively. In Fig. [2] we fix Q 2 at three different values: Q 2 = 0.3 GeV 2 in 
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FIG. 2: Nonsinglet SF F 2 (x, Q 2 ) vs x at (a) LO and (b) NLO. The energy scale is Q 2 = 0.3 GeV 2 in (a.l), 
(b.l), Q 2 = 1 GeV 2 in (a.2), (b.2) and Q 2 = 100 GeV 2 in (a.3), (b.3). The solid line represents the FAPT 
results and the dashed line - the pQCD ones. 


(a.l), (b.l), with the initial point Q 2 = 1 GeV 2 in (a.2), (b.2), and Q 2 = 100 GeV 2 in (a.3), 
(b.3), where again (a) and (b) represent the LO and NLO, respectively. 


V. Summary 


The main goal of this work is to propose a new theoretical tool for the DIS analysis, based 
on Fractional Analytic Perturbation Theory, to the DIS community. This approach allows 
one to analyze formally the leading-twist structure function in the whole Q 2 range. This 
conclusion is explicitly shown in Figs. |T]and[2j The scheme of the approach is formulated 
in Sec. m and applied for data processing in Sec. [TV] Our consideration is restricted to the 
leading twist. The higher twist contributions (HT) can be taken into account by a fit of 
experimental data together with PDFs. Moreover, the role of the stability of APT for this 
fit was pointed out in [ll], 27-2f| (and in Introduction here). Our investigation reveals the 
following main features of applying FAPT: 


• Structure function F 2 (x, Q 2 ) at fixed x changes very slowly in the entire range of Q 2 . 

• At high Q 2 evolution ( Q 2 > 100 GeV 2 ) the pQCD and FAPT distributions become 
practically equal. 


• The evolution in FAPT is more gradual (i.e., it evolves slower) and smoother than in 

pQCD. 


• The new analytic (FAPT) series converge faster than the pQCD series. From inspec¬ 
tion of Figs. [HE] it is obvious that the one- and the two-loop FAPT approximations 
do not differ significantly from each other (the difference is less than 1%). 


In this work, we have analyzed only the nonsinglet part, the consideration of the singlet 
part can be performed along the same line but requires more complicated formulas and 
cumbersome numerical calculations. This is the task for forthcoming investigation. Other 
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important issues to complete this FAPT approach as the reliable tool for DIS is to add the 
target mass corrections (TMC) and the aforementioned HT contributions in our scheme of 
calculation. These improvements will help one to clarify in future the behavior at very low 
energies ( Q 2 ~ 0.3 GeV 2 ) in more detail, ft would be important to emphasize, that the 
FAPT approach admits investigation of the HT contributions in the most sensitive regime 
of moderate/small Q 2 due to the high stability of the radiative corrections. 
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Appendix A: Explicit expressions for NLO /3-functions, anomalous dimensions and 
coefficient functions of DIS (nonsinglet case) 


The renormalization group equation for a s 
up to the NLO approximation is given by 


a s (L) 

47T 


at the expansion of the /3-function 


d_ 

dL 


a s (L) 


-PML)) 


where the first two beta coefficients are 


—/3 0 a 2 (L) — /3i a 3 s (L) + ... , 


(Al) 


11 4 34 / 20 \ 

A) = y C A ~ ~ T^Jif , /A = — C\ - F 4Cf + — Ca ) T R n f . (A2) 

The anomalous dimensions of composite operators in LO, 7 ns( n ), NLO 7 ns( n ) an d the 

coefficient function C’/jg (n) are expressed by means of transcendental sums S a (n ), see, e.g., 

H], 



2 C F 



2 

n{n + 1) 


+ 4 (Si(n) — 1) 


(A3a) 
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(V^c a )x{i6S iW A^±D + 1( 

S 2 (n) - S± (^j + 645 ± (n) + 24 S 2 (n) - 3 - 8 Sf 

Q O ^ O . ^ 


n(n + 1) 


, 3n 3 + n 2 — 1 2n 2 + 2n + 1 

T 6 - 3 (n + l) 3 

536 2n + 1 

+ ■" 


n 3 (n + 1) J n 


+C F C A { S\ (n) 
+S 2 (n 


9 n 2 (n + l) 2 


_52 8 

3 n(n + 1) 6 


H 4 


- 16S'i(n)S' 2 (n) 

151n 4 + 263n 3 + 97n 2 + 3n + 9 
9n 3 (n + l) 3 
11 n 2 + 5n — 3 


,ri ,T rr, f 160 c ( , , 32 c ^ , 4 , ic lln 2 + 5n- 
+C F N F T R S^n) + -S 2 (n) + - + 16 ^ + J)2 


(A3b) 


4s (n) = C F ( 25' 2 (n) + 3 S^n) - 2 S 2 (n) - 


2Si(n) 
n(n + 1) ' n ' n + 1 


3 

+ - + 


+ ^-9 ,(A3c) 


n z 


On the other hand, the series S a (n ) = ^fc=i can be expressed via the generalized Riemann 
£ functions, see 391], that are analytic functions in both variables a,n: 


Si(n) = 0(n + 1) - 0(1), (A4a) 

S 2 (n) = C(2) — 4>'(n + 1) = C(2) — C(2, n + 1), (A4b) 

S a {n) = ((a) -((a,n + 1). (A4c) 


r OO / - \ ^. 

^(n) = S- 2,1 = —gC(3) T (b + n) 2 ^ + n + ^ _ V’C 1 )) ■ 

fc=i ^ + ' 


(A4d) 


For the S ± and S series we use the notation given in |40|, |4lj. Performing the analytic 
continuation from even n, S£, and from odd n, S~ ( see for details 0) one obtains 


S+(n/ 2) ->2“- 1 [S a (n) + S+ a (n)] =2“" 1 [S a (n) + C(«) - 4>(-l, a, n + 1)] - C(«), (A5a) 
S~(n/ 2)^2 Q ~ 1 [S Q (n) + 5l a (n)]=2“- 1 [5 a (n) + C(«) + *(-l,a,n + 1)] - C(«), (A5b) 

where <&(z,a,v) is the Lerch transcendent function [39]. The expressions on the r.h.s. of 
Eqs. (IA5p are now analytic functions in both variables a,n - this is a new result. 


Appendix B: QCD evolution of moments up to NLO 

The PDFs are the nonsinglet fas{x, Q 2 ) an d singlet fs(x,Q 2 ) parton distribution func¬ 
tions, 

hs(x,Q 2 ) = u v (x,Q 2 ) - d v (x,Q 2 ), (Bl) 

fs(x, Q 2 ) = u v (x, Q 2 ) + d v (x, Q 2 ) + S{x, Q 2 ) = V(x, Q 2 ) + S(x, Q 2 ), (B2) 

whereas V(x, Q 2 ) is the distribution of valence quarks and S(x,Q 2 ) is the sea quark distri¬ 
bution. More generally, the NS PDF is a combination of the forms u — d and d — u but for 
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our consideration we focus on the nucleon scattering provided by combination (Bl). The 
moments representation for PDFs is defined as 

hs(n,Q 2 ) = [ dxx n ~ l fas(x,Q 2 ), (B3) 

Jo 

fs(n,Q 2 ) = [ dxx^fsix^ 2 ). (B4) 

Jo 

The moments Mns(?L y 2 ) for the structure function F^s{x,y 2 ) follow from the Mellin con¬ 
volution F d = C * f, 


Fns(z, h 2 ) 
A/ns(?L h 2 ) 


(Cns * /ns) (z,fi 2 ) = / C N s(2/,a s )/ NS (>,/i 2 ) 5{z - x 
Jo 

C'ns(^) a s (y 2 )) ‘ /ns Ob y 2 ) • 


y)dydx, (B5a) 
(B5b) 


Here Cns <b) is the nonsinglet coefficient function of the process that can be presented as 
the perturbation series C^s(x,a s ) = 1 + «. s (Q 2 )Cxs(x) + 0(a 2 ); ( n ) i n Appendix lAl is 


the moment of the C'^g (x). The QCD evolution of the moments Mns up to NLO of is given 
by (see Ref [38]) 


ATns (n, Q 2 ) — 


1 + ^ns (n)a s (Q 2 ) /l + (fii/j3o)a s (Q 2 ) 


2 \ \ p(n) r 


1 + <?nsVH(<2o) V 1 + (PMhs(Qo) 


where 

and: 


dm(n) = 7Ns( n )/2A)> p(n) = 


„(!)/ 


„(o). 


A 


Po 


a s (Q 2 ) 

d NS (n) 


a s (Q 2 0 )_ 



xM n s(n, Ql), 

(B6a) 

), 


(B6b) 

n)\ 

/ 


(B6c) 


The coefficients of anomalous dimension in LO and NLO and the coefficient function in NLO 
are given in Eqs. flA3a[) . Appendix [Al In the case of the nucleon structure function F 2 (x, Q 2 ), 
one needs to take into account only even values of n in the NLO anomalous dimension. 


Appendix C: Accuracy of the Jacobi Polynomial method 


The accuracy of the evaluation of the structure functions depends on the method we use; 
therefore, it is indispensable to verify it in our approach. The Jacobi Polynomial method 
promises us a good enough accuracy for the evolution, as was shown in previous works 

H). 


see 


This method is applied directly to the terms of the Bjorken variable x, but it affects the 
Q 2 -dependence indirectly. Therefore, it is necessary to confirm the x-range applicability of 
the JP method. To this end, we compare the results of the JP approach with the “exact” 
numerical calculations of inverse Mellin moments following Eq. (12. 6 p but only in the one-loop 
approximation due to technical limitations. The comparison of these two results in Fig. [3] 
demonstrates a very good accuracy. So, in order to clarify it, we perform a zoom in x, going 
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to a lower x-region (~ 10 -2 ). We see in Fig. [4] that the JP method gradually loses precision 
starting at x < 0.02. Also, we can see that in this range, the difference between these two 
methods reaches 5%. 


Appendix D: Accuracy of the rational approximation 

Here we investigate the accuracy of the rational approximation for the two-loop evolution 
factor 


m(n, Q 2 ) 


f i + (A/A) ) Q s(Q 2 ) \ p( ' n ' > 

\1 + {Pi/Po) a s{Qo) / 


(Dl) 


in Eq. (l3.12j) for pQCD, and Eq. (l3.13p for FAPT, respectively. The expansion of the factor 
in power series up to NLO leads 


™pQCd (™,Q 2 ) - 

m pQCD{n,Q 2 ) - 


1 + (A/^ 0 )p(w)a s (Q 2 ) 

1 + (Pi/Pq )p{n)a s {QlY 

1 + (Pi/Po)p(n)a s {Q 2 ) + {Pl/2Pl)p{n){p{n) - 1 )a 2 s (Q 2 ) 

1 + {Pi/Po)p(n)a a {QD + {131/2Pl)p{n){p{n) - l)a2(Q 2 )’ 


(D2a) 

(D2b) 


where m^{n,Q 2 ) and 7n^ 2 \n,Q 2 ) represent the approximation up to 0(a s ) and 0(a 2 ), 
respectively. The corresponding “FAPT form” of (ID2D is given by 


^FAPT^b Q ) — 

m FAPT ( n iQ") — 


1 + (Pi/ f3o)p(n)-Ai{Q 2 ) 

1 + Wi/Po)p{n)A 1 {QlY 

1 + WMpHA^Q 2 ) + (/ 3j/2/3$)p{n){p{n ) - 1 )A 2 {Q 2 ) 

1 + (Pi/Po)p(n)Ai(Qo) + (Pi/2Pi)p{n){p{n) - l)A 2 {Ql)' 


(D3a) 

(D3b) 


Combing the approximations in Eqs. (ID21 ID3|) in quantity A m^ 12 ' 1 = \m/ 1 ^ — m^\/m^ we 
obtain the accuracy better than 1% for any n < 13 (since JP expansion contains only 13 
terms for good approximation), in both cases of pQCD and FAPT for two different ranges 


0.10 
0.08 
0.06 
0.04 
0.02 

0.0 0.2 0.4 0.6 0.8 1.0 °' 0 °0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 3: Nonsinglet SF F 2 {x) vs x in LO (a), at the energy scale Q 2 = 0.3 GeV 2 and (b) Q 2 = 100 GeV 2 . 
The solid (red) line represents the FAPT result and the dashed (blue) line the pQCD one in the JP method. 
The thick squares (red) and spheres (blue) represent the result of the “exact” numerical Inverse Mellin 
transform. 
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FIG. 4: Nonsinglet SF F^x) vs x at LO , at energy scale (a) Q 2 = 0.3 GeV 2 , and (b) Q 2 = 100 GeV 2 . 
The solid (red) line represents the FAPT outcome and the dashed (blue) line the pQCD one in the JP 
method. The thick squares (red) and spheres (blue) represent the result of the “exact” numerical Inverse 
Mellin transform. 


TABLE I: The accuracy in per cent of the difference of the approximations Am^ 12 ' = 

If 

pQCD: AmpQCD and f° r FAPT: Am^p T . The results are presented in two ranges of Q 2 : low Q 
and high Q 2 ~ 100 GeV 2 . 


I rrS 1 ^ — m^ I 


(i) 


for 


1 GeV 2 3 4 5 6 


n 

2 

4 

6 

8 

10 

12 

^ m pQCD% 

0 . 74 / 

0 . 49 / 

0.26 / 

0 . 06 / 

0 . 12 / 

0.3 

Q 2 ~ 1 GeV 2 /100 GeV 2 

0.65 

0.45 

0.25 

0.06 

0.12 

0.3 

Am,p^p T % 

0.03 / 

0.02 / 

0.01 / 

0.00 / 

0.01 / 

0.01 

Q 2 ~ 1 GeV 2 /100 GeV 2 

0.18 

0.13 

0.07 

0.02 

0.04 

0.09 


of energy, i.e., low Q 2 ~ 1 GeV 2 and high Q 2 ~ 100 GeV 2 . The results collected in Table 
HI demonstrate that for both ranges of energy FAPT has a better convergence than pQCD; 
even more, the accuracy is improved for FAPT at low Q 2 r\j 1 GeV 2 (really, pQCD must 
be worse but we have a low starting point Qq — 1 GeV 2 ). The strong hierarchy of FAPT 
couplings, |^ F j A L PT) (Q 2 )| -C |^ FAPT) (Q 2 )|, remains valid even at very low \Q 2 \, cf. [loj]. 


[1] D. V. Shirkov and I. L. Solovtsov, |hep-ph/9604363|; Phys. Rev. Lett. 79, 1209 (1997) 
j[hep-ph /9704333]. 

[2] K. A. Milton, I. L. Solovtsov Phys. Rev. D 55, 5295 (1997) hep-ph/9611438]. 

[3] K. A. Milton, I. L. Solovtsov and O. P. Solovtsova, Phys. Lett. B 415, 104 (1997) 
arXiv:hep-ph/97064091. 

[4] K. A. Milton, O. P. Solovtsova Phys. Rev. D 57, 5402 (1998) |hep-ph/9710316|. 

[5] D. V. Shirkov, Theor. Math. Phys. 127 , 409 (2001) jhep-ph/0012283]; 119 , 438 (1999). 

[6] A. I. Karanikas and N. G. Stefanis, Phys. Lett. B 504, 225 (2001) [Erratum-ibid. B 636, 330 
(2006)] [hep-ph/0101031]]. 


14 















[7] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Rev. D 72, 074014 (2005) [Erratum- 
ibid. D 72, 119908 (2005)] [hep-ph/0506311]. 

[8] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Rev. D 75, 056005 (2007) [Erratum- 
ibid. D 77, 079901 (2008)] (hep-ph/0607040j. 

[9] A. P. Bakulev, S. V. Mikhailov, N. G. Stefanis, JHEP 1006, 085 (2010) [arXiv: 1004.4125 ]. 

[10] A. P. Bakulev, Phys. Part. Nucl. 40, 715 (2009) [arXiv:0805.0829 ] (arXiv preprint in Russian); 
N. G. Stefanis, Phys. Part. Nucl. 44, 494 (2013) [arXiv:0902.4805]. 

[11] A. V. Sidorov and O. P. Solovtsova, Mod. Phys. Lett. A 29 (2014) 1450194, [arXiv: 1407.6858]; 
Nonlin. Phenom. Complex Syst. 16, 397 (2013), [arXiv:1312.3082 ]. 

[12] D. V. Shirkov and I. L. Solovtsov, Theor. Math. Phys. 150, 132 (2007) |hep-ph/0611229|. 

[13] A. V. Nesterenko, Phys. Rev. D 62, 094028 (2000); Phys. Rev. D 64, 116009 (2001); Int. J. 
Mod. Phys. A 18, 5475 (2003). 

[14] A. V. Nesterenko and J. Papavassiliou, Phys. Rev. D 71, 016009 (2005); A. C. Aguilar, 
A. V. Nesterenko and J. Papavassiliou, J. Phys. G 31, 997 (2005). J. Phys. G 32, 1025 (2006) 
|hep-ph/0511215|; A. V. Nesterenko, arXiv:0710.5878, 

[15] A. I. Alekseev, Few Body Syst. 40, 57 (2006) |hep-ph/0503242j|. 

[16] Y. Srivastava, S. Pacetti, G. Pancheri and A. Widorn, In the Proceedings of e + e~ Physics 
at Intermediate Energies, SLAC, Stanford, CA, USA, 30 April - 2 May 2001, pp T19 
[hep-ph/0106005]. 

[17] B. R. Webber, JHEP 9810, 012 (1998) jhep-ph/9805484]. 

[18] G. Cvetic and C. Valenzuela, J. Phys. G 32, L27 (2006) |hep-ph/0601050|. 

[19] G. Cvetic and C. Valenzuela, Phys. Rev. D 74, 114030 (2006) |hep-ph/0608256[. 

[20] G. M. Prosperi, M. Raciti and C. Simolo, Prog. Part. Nucl. Phys. 58, 387 (2007) 
[hep-ph/06072091. 

[21] G. Cvetic and C. Valenzuela, Braz. J. Phys. 38, 371 (2008) [arXiv:0804.0872 ]. 

[22] Y. O. Belyakova and A. V. Nesterenko, Int. J. Mod. Phys. A 26, 981 (2011) [arXiv:1011.1148j. 

[23] G. Cvetic, R. Kogerler and C. Valenzuela, J. Phys. G 37, 075001 (2010) [arXiv:0912.2466j. 

[24] G. Cvetic, R. Kogerler and C. Valenzuela, Phys. Rev. D 82, 114004 (2010) [arXiv: 1006.4199]. 

[25] C. Contreras, G. Cvetic, R. Kogerler, P. Kroger and O. Orellana, arXiv:1405.5815 . 

[26] C. Ayala, C. Contreras and G. Cvetic, Phys. Rev. D 85, 114043 (2012) [arXiv: 1203.6897 ]. 

[27] A. V. Kotikov, V. G. Krivokhizhin and B. G. Shaikhatdenov, Phys. Atom. Nucl. 75, 507 (2012) 
[arXiv: 1008.0545 [hep-ph]]; G. Cvetic, A. Y. Illarionov, B. A. Kniehl and A. V. Kotikov, Phys. 
Lett. B 679, 350 (2009) [arXiv:0906.1925 [hep-ph]]. 

[28] R. S. Pasechnik, D. V. Shirkov and O. V. Teryaev, Phys. Rev. D 78, 071902 (2008) 
[arXiv:0808.0066 [hep-ph]]; R. S. Pasechnik, D. V. Shirkov, O. V. Teryaev, O. P. Solovtsova and 
V. L. Khandramai, Phys. Rev. D 81, 016010 (2010) [arXiv:0911.3297 [hep-ph]]; V. L. Khan- 
dramai, R. S. Pasechnik, D. V. Shirkov, O. P. Solovtsova and O. V. Teryaev, Phys. Lett. B 
706, 340 (2012) [arXiv:1106.6352], 

[29] O. Teryaev, Nucl.Phys.Proc.Suppl. 245 (2013) 195 [arXiv: 1309.1985 ]. 

[30] P. Allendes, C. Ayala and G. Cvetic, Phys. Rev. D 89, 054016 (2014) [arXiv:1401.1192j. 

[31] V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972) [Yad. Fiz. 15, 781 (1972)]; 
L. N. Lipatov, Sov. J. Nucl. Phys. 20, 94 (1975) [Yad. Fiz. 20, 181 (1974)]; G. Altarelli and 
G. Parisi, Nucl. Phys. B 126, 298 (1977); Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977) 
[Zh. Eksp. Teor. Fiz. 73, 1216 (1977)]. 

[32] V. G. Krivokhizhin et ah, Z. Phys. C 36, 51 (1987); V. G. Krivokhizhin et ah, Z. Phys. C 48, 
347 (1990). 


15 


[33] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63, 189 (2009) 
[arXiv:0901.0002 ]. 

[34] A. J. Buras, Rev. Mod. Phys. 52, 199 (1980). 

[35] V. G. Krivokhizhin and A. V. Kotikov, Phys. Atom. Nucl. 68 (2005) 1873 [Yad. Fiz. 68 (2005) 
1935]; V. G. Krivokhizhin and A. V. Kotikov, Phys. Part. Nucl. 40, 1059 (2009). 

[36] A. P. Bakulev and V. L. Khandramai, Comput.Phys.Commun. 184(2013) 1, 183 

arXiv: 1204.2679. 

[37] C. Ayala and G. Cvetic, Comput. Phys. Commun. 190, 182 (2015) [arXiv: 1408.6868 [hep-ph]]; 
arXiv:1411.1581 [hep-ph], 

[38] F. J. Yndurain, The Theory of Quarks and Gluons Interactions (Fourth Edition) (Springer- 
Verlag, Berlin, 2006). 

[39] A. Erdelyi, W. Magnus, F. Oberhettinger and F. G. Tricomi (1953), Higher Transcendental 
Functions. Vol. I, McGraw-Hill Book Company, Inc., New York-Toronto-London. 

[40] J. Blumlein and S. Kurth, Phys. Rev. D 60, 014018 (1999); 

[41] A.V. Kotikov and V. N. Velizhanin, hep-ph/0501274 

E. G. Floratos, D. A. Ross and C. T. Sachrajda, Nucl. Phys. B 129, 66 (1977) [Erratum-ibid. 
B 139, 545 (1978)]. A. Gonzalez-Arroyo, C. Lopez and F. J. Yndurain, Nucl. Phys. B 159, 
512 (1979). G. Curd, W. Furmanski and R. Petronzio, Nucl. Phys. B 175, 27 (1980). 


16 



